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FOIL  FOCUSING  OF  ELECTRON  BEAMS 


I.  INTRODUCTION 

Accelerating  intense,  high-current  electron  beams  to  ultrarelativistic 

energies  requires  a  focusing  mechanism  to  overcome  the  defocusing  from  beam 

emittance  and  space-charge.  The  need  for  focusing  is  most  critical  in  the 

initial  accelerator  stages  when  the  beam  energy  y  is  smallest.  As  y  rises, 

the  need  for  focusing  decreases  because  the  azimuthal  magnetic  field  of  the 

_2 

beam  offsets  the  defocusing  space-charge  field  to  order  y  ,  and  because 

-1  .  1 

the  unnormalized  beam  emittance  decreases  as  y  in  an  ideal  accelerator. 

The  two  main  focusing  methods  used  in  high-current  linear  induction 

accelerators  are  solenoidal  fields  and,  more  recently,  radial  electrostatic 

2 

fields  from  an  ion  column.  Solenoidal  focusing  requires  precise  magnetic 
alignment  and  a  large  investment  in  magnetic  energy.  Ion  focusing,  on  the 
other  hand,  requires  a  plasma  column  with  the  attendant  problems  of  gas 
pumping,  plasma  instabilities,  loading  of  the  accelerating  gaps,  etc.  Both 

transport  techniques  are,  moreover,  subject  to  beam  instabilities,  either 

3  4 

from  cavity  modes  or  from  ion  motion. 

5  , 

An  alternative  transport  technique,  first  suggested  by  Adler,  is  a 
series  of  thin  conducting  foils  placed  transverse  to  the  beam  axis.  Foils 
act  like  thin  electrostatic  lenses  that  focus  the  beam  via  image  charges. 

Foils  are  easy  to  implement  and  may  help  stabilize  the  beam.^’^  Initial 

6  8  9 

experimental  studies  of  foil  focusing  were  encouraging,  although  Meger 

has  reported  difficulty  at  beam  currents  above  10  kA. 

In  this  paper  we  model  foil  focusing  using  the  thin-lens  approximation 
5  10 

assumed  by  Adler  and  Humphries.  With  this  model  we  determine  the 
dependence  of  foil  focusing  on  the  beam  parameters,  study  its  effect  on 
beam  emittance,  and  examine  stability  of  multi-foil  transport.  Although 
the  thin-lens  approximation  is  generally  valid  only  for  paraxial  beams  of 
current  1^  «  17y  kA,  preliminary  numerical  simulations  suggest  that  the 
analytical  predictions  remain  roughly  intact  even  when  the  paraxial 
approximation  fails.  The  numerical  simulations  will  be  presented  in  a 
later  publication. 
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The  .paper  .is  organized  as  follows.  The  focusing  properties  of  an 
isolated  thin  foil  are  determined  in  Sec.  II  for  several  different, 
axisymmetric  beam  profiles.  Emittance  growth  from  anharmonic  foil  focusing 
is  computed  in  Sec.  III.  Multi-foil  transport,  including  beam  stability 
and  electrical  interactions  between  adjacent  foils,  is  considered  in  Sec. 
IV.  Apertures,  foil  scattering,  plasma  effects,  and  beam  centering  are 
covered  briefly  in  Sec.  V.  Section  VI  summarizes  the  results. 

II.  FOIL  FOCUSING 

A.  Electrostatic  Potential  from  a  Single,  Thin  Foil 

A  conducting  foil  placed  perpendicular  to  the  path  of  an  electron  beam 
induces  foil  currents  and  charges  that  short  the  electromagnetic  fields  of 
the  beam.  The  foil  charges  produce  electrostatic  fields  that  persist  well 
outside  the  foil  and  act  to  focus  the  beam.  The  foil  currents,  on  the 
other  hand,  little  affect  the  beam  because  they  are  closely  spaced  and  flow 
anti-parallel  on  the  two  faces  of  the  foil.  The  primary  effect  of  the 
foil,  therefore,  is  to  focus  the  beam  electrostatically.'*''*'  In  this  section 
we  compute  the  electrostatic  foil  fields  for  several  different  beam 
profiles. 

To  model  foil  focusing,  consider  an  ultrarelativistic,  paraxial  beam 
propagating  in  a  conducting  pipe  of  radius  b.  The  beam  is  paraxial  in  that 
the  beam  electrons  follow  nearly  straight  and  parallel  orbits  over  axial 
distances  ~  b.  This  assumption  allows  us  to  treat  the  foil  as  a  thin  but 
imperfect  lens.  In  addition,  we  neglect  foil  scattering  and  assume  that 
neither  the  beam  current  nor  radius  changes  rapidly  within  a  time  b/c.  The 
last  assumption  allows  us  to  neglect  inductive  and  time-retardation 
corrections  to  the  fields. 

The  foil  charge  produces  an  electrostatic  potential  satisfying 

=  0,  (1) 

with  <}>£  equal  to  zero  at  z  =  ±°>  and  on  the  pipe,  r  =  b.  On  the  foil 
itself, 

♦f  =  -  V  (2) 
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where  $  is  the  potential  from  the  beam  (and  associated  plasma,  if  any)  in 
the  Absence  of  the  foil.  The  foil  thus  converts  the  ambient  transverse 
field,  E^o  =  *nt0  a  longitudinal  field,  Ez^  =  -3<t>£/3z.  The 

longitudinal  foil  field  E^  accelerates  the  beam  as  it  approaches  the  foil 
and  decelerates  the  beam  as  it  leaves  the  foil.  The  major  influence  on 
beam  dynamics,,  however,  is  the  electrostatic  pinch  caused  by  The 

longitudinal  field  Ez£  is  much  less  important  because  it  is  anti-symmetric 
about  the  foil,  and  because  the  transverse  energy  of  a  paraxial  beam  is 
much  smaller  than  the  longitudinal  energy. 

For  an  axisymmetric  beam,  a  single  foil  located  at  z  =  0  generates  an 
electrostatic  potential  that  can  be  expressed  as 

00 

<f>f(r,z)  =  I  An  expC-xJzJ/b)  J^^r/b),  (3) 

n=l 

where  ^  is  the  nth  root  of  the  Bessel  function  JQ  and 

b 

Jdr  r  JQ(Xnr/b)  $o(r,0)  (4a) 

0 

b 

v  h  ,‘fy  r  Jl<XnC/b>  Ero(r’°>  C4b) 

V  JrV  o 

b 

■  -  70777  Jdr  r  Jo<Xnr/b>  &  If  r  Ero(r'0)]-  <4c> 

*n  W  0 

Here  E  =  -3<J>  /3r  is  the  radial  electrostatic  field  from  the  beam  and 
ro  0 

plasma,  and  we  have  integrated  by  parts  using  the  Bessel  identities, 

Jx(x)  =  -  dJQ/dx  and  xJQ(x)  =  d(xJ^)/dx. 

B.  Radial  Impulse  Produced  by  a  Foil 

As  Eq.  (3)  demonstrates,  the  foil  pinch  field,  E  -  =  -Sf^/Sr,  is  strong 

only  within  an  axial  distance  ~  b/X^.  In  this  short  distance,  the  foil 

alt*  .s  the  radial  momentum  but  not  the  orbital  radius  of  the  paraxial 

electrons.  Assuming  no  change  in  radius  r  or  axial  velocity  v  ,  the 

z 

electrons  acquire  a  radial  impulse  given  by 


An  "  “  ,2 .2 


b  Ji<V 
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Here  we  have  set  vz  equal  to  the  speed  of  light,  c.  The  assumption  of 
constant  r  and  constant  v constitutes  the  thin-lens  approximation. 

To  quantify  the  analysis,  consider  a  rigid-rod  beam  and  a  plasma  that 
neutralizes  the  beam  charge  by  a  fraction  f£(r)  <  1.  Most  applications 
operate  in  vacuum  where  f  =  0,  but  for  generality  we  allow  f  t-  0.  For  a 
long  rigid-rod  beam,  3/3z  =  0  and  the  ambient  space-charge  field  from  the 
beam  and  plasma  satisfies 


1  3_ 
r  3r 


4n 

c 


<i-fc>V 


(6) 


where  r )  is  the  beam  current  density.  Inserting  Eq.  (6)  into  Eq.  (4c) 
yields 


We  now  compute  An  for  four  different  beam  profiles,  assuming  constant  f£. 

The  influence  of  the  foil  on  f  is  discussed  later. 

c 

For  a  flat-topped  beam  with  constant  out  to  a  radius  r^  <  b,  Eq.  (7) 
reduces  to 


<1-fc>Ibb  Jl(Xnrb/b> 
c  rb  *n3 


(3a) 


2 

where  1^  =  itr^  is  the  beam  current.  Inserting  Eq.  (8a)  into  the  impulse 
equation  (5)  yields  Adler' sb  result  for  the  deflection  angle  from  the  foil: 


S0f  = 


,,  b.  r  ¥Vb/b>  Jl<V/b> 

■  JA  rb  ll  ^  jfo) 


3 

where  P2  =  yinc  is  the  axial  momentum,  IA  =  yinc  /e  =  17y  kA  is  the  Alfven 

current,  and  y  is  the  Lorentz  factor.  The  mean-squared  radius  for  a  flat- 
2  2 

topped  beam  is  R  =  r,  /2. 

D  2  2 
For  a  Gaussian  beam  with  varying  as  exp(-r  /R  ), 
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f  2  )  (1"fc):tb  2 

An  =  -  fcvvJ  - c  expt-O^R/Zb)2], 


(1-f  )I.  •  expl-(X_R/2b)/] 

S6  =  -  8  - - 5— | - J1(Xnr/b>' 

n=l  X„2  J^) 


to  order  exp(-b*7R  ).  The  rms  radius  equals  R. 


2,_  2, 


For  a  parabolic  beam  with  varying  as  (1-r  /r^  )  for  r  <  r^, 


An  =  - 


<1-fc>Ib  f  8b/rb 


8b/r,  r  f  . 

^Tw\ 


Jo(Xnrb/b) 


,  (10a) 


50f  =  - 


(1-f  )I,  ®  f  8b/r,  )*  . 

r-^I  ~2 - —  Jl<V/b)  FT  Jl(Xnrb/b) 

A  n=l  IXjj  J^X^J  L  nCb 


Jo(Xnrb/b) 


2  2 

The  mean-squared  radius  is  R  =  r,  /3« 


An  =  -4 


v-Vh  fxif  VWb> 


l*nJ  J2^)  [Xj2  -  (X^/b)2]’ 


(10b) 


For  a  Bessel  profile  with  varying  as  JQ(X]Lir/rtJ)  for  r  <  » 


(11a) 


„  „  ;  fk)  J,<Wb>  ¥V/b> 

*  h  n=lW  J12(X„)  IXj2  -  (XnL'b/b)2j 

2  2  2  2 
The  mean-squared  radius  is  R  =  (1-4/X^  )r^  -  (0.555  r^)  . 


(lib) 


C.  Foil  Focal  Length 

In  Fig.  1  we  plot  the  foil  focal  length, 


1 


(12) 


ft(r)  *  -  r/S9f(r), 

for  three  values  of  b/R  for  each  of  the  four  profiles.  For  all  profiles, 
f^  is  proportional  to  YR/(l-Jtc)I^>  and  it  increases  with  particle  radius  r 
but  decreases  with  wall  radius  b.  The  dependence  on  r  means  that  the  lens 
is  imperfect,  with  the  beam  image  distorted  and  emittance  degraded.  Any 
spread  in  particle  energy  in  a  given  transverse  segment  of  the  beam  further 
degrades  the  emittance  through  the  dependence  of  f^  on  y*  The  variation  in 
f^  with  r  is  the  least  for  a  flat-topped  beam,  for  which  the  self-field  ErQ 
is  harmonic  (proportional  to  r  for  r  <  r^),  and  is  the  greatest  for  the 
Gaussian  profile  which  extends  to  large  r.  The  parabolic  and  Bessel 
profiles  both  cut  off  sharply  and  produce  nearly  identical  f^. 

The  dependence  of  the  focal  length  on  the  beam  parameters  makes  foil 
focusing  difficult  to  control,  but  it  also  provides  a  means  for  altering 
the  beam  envelope,  R(C)1;  here  C  =  ct-z  measures  distance  into  the  beam 
pulse.  Controlled  manipulation  of  R(0  is  possible  only  in  regions  where 
I^/y  varies  in  a  known  manner  with  t,.  Unpredictable  fluctuations  in  1^,  y> 
or  R  produce  like  fluctuations  in  f^  that  can  grow  with  propagation 
distance  z.  See  Sec.  IV  for  further  discussion. 

The  paraxial  approximation  is  valid  only  if  f^  »  b,  or  equivalently, 
S0j  «  1.  For  foils,  this  requires  that  (l-f^I^  «  IA>  At  high  beam 
currents,  1^  >  0.3  I^/(l-fc),  foil  focusing  becomes  so  strong  that  the 
axial  and  transverse  particle  velocities  become  comparable.  Beam  quality 
is  then  poor,  although  efficient  foil  transport  is  possible  even  at 
currents  above  the  space-charge  limiting  current. ^ 

III.  BRAM  EMITTANCE 


A.  Imperfect  Lenses 

A  useful  parameter  for  characterizing  beam  quality  is  the  rms 
transverse  emittance.  For  simplicity,  we  consider  an  axisymmetric. 

12 

non-rotating,  monoenergetic  beam  with  a  normalized  emittance  defined  by 

e  =  Y  [R2<ex2>  -  <r1.gi>2]1/2,  (13) 
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where  the  mean-squared  beam  radius  is  given  by 
* 

R2  =  <rx2>,  '  (14) 

and  r^  and  g^  =  dr^/dz  are  the  transverse  position  and  velocity, 
respectively,  of  a  beam  electron.  The  angular  brackets  denote  an  average 
taken  over  all  particles  in  a  given  transverse  segment.  Harmonic  forces 
preserve  e,  but  anharmonic  forces  generally  do  not.*  In  an  ideal 
accelerator  with  harmonic  forces,  the  unnormalized  emittance  e/y  determines 
beam  expansion  and  decreases  as  y  rises. 

To  compute  the  emittance  change  produced  by  an  anharmonic  thin  lens,  we 
express  the  transverse  velocity  of  each  beam  electron  as 

gA  =  grr  +  (15) 

where  Sg^  is  a  random  thermal  velocity  and  Pr  is  the  radial  fluid  velocity 
defined  by 

Pr(r'C)  =  JJd2§JL  §|/*  (16) 

Here  F(r^,g^,C)  is  the  electron  distribution  function,  and  r  is  the  unit 
radial  vector.  The  thermal  component  satisfies  by  definition 

Jjd\  Sgj L-r  F(ri,gi,C)  =  0.  (17) 

A  thin  lens  alters  g^  by  60^  r: 

gA  *  g±  +  S0f  r  =  Sg^  +  [gr  -  (r/f£)]r.  (18) 

According  to  Eqs.  (13)-(17),  this  alters  the  square  of  the  emittance  by 

S(e2)  =  ef2  +  ef2,  (19) 

where 

efo  5  Y*  tR2«r/f())2>  -  <r2/fe>2] ,  (20a) 
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(20b) 


and 

ef2  *  2 r2  l<r2/ft>  <Prr>  -  R2 
2 

The  emittance  change  is  independent  of  the  beam  velocity  distribution 
and  is  always , positive;  by  contrast,  e^2  is  controlled  by  changes  in  the 
beam  profile  and  can  be  positive  or  negative.  In  deriving  Eq.  (19),  we 
have  used  the  fact  that 


<g(r)  5g±*r>  =  0,  (21) 

for  all  g(r). 

Several  points  are  worth  noting  from  results  (19)  and  (20).  First,  a 

perfect  lens  does  not  alter  the  beam  emittance:  =  5e  =  0  for 

3f^/9r  =  0.  Second,  the  change  in  emittance  from  an  imperfect  lens  is 

independent  of  the  beam  thermal  velocity  a  hot  beam  thus  suffers  the 

same  emittance  change  as  a  cold  beam,  assuming  identical  fluid  velocity 

12 

Pr(r>0*  Third,  if  the  beam  evolves  self-similarly,  the  term  equals 

zero  and  the  emittance  squared  increases  by  ef2  >  0;  self-similar  evolution 

ro  2 

means  that  the  beam  profile  is  invariant,  with  ur  «  r/R  and  J^R  a  function 
of  r/R  only.  An  anharmonic  thin  lens  thus  increases  the  square  of  the 
emittance  by  e^,  unless  the  beam  changes  its  profile  vhile  striking  the 
lens.  We  defer  discussion  of  profile  changes  until  Sec.  III-C. 


B.  Scaling  of  e^o 

The  emittance  increase  can  be  estimated  using  a  generic  model  for 
the  foil  focal  length: 


f£(r)  ~  goR  (l-£c)Ib 


(1  +  gjr/R2), 


(22a) 


where  gQ  characterizes  the  lens  strength  and  g^  is  an  anharmonic  lens 
coefficient.  This  form  fits  the  four  profiles  in  Fig.  1  well  for  r  <  2R. 
Table  1  lists  the  approximate  coefficients  gQ  and  g^  for  the  flat-topped, 
Gaussian,  parabolic,  and  Bessel  profiles  for  several  values  of  b/R.  The 
anharmonic  coefficient  g^  is  smallest  for  the  flat-topped  profile  and 
largest  for  the  Gaussian  profile. 
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TABLE  1:  VALUES  OF  (gQ, 


Beam  Profile 

b/R  =  2 

5 

10 

Flat-topped 

(0.45,  0.12) 

(0.37,  0.13) 

(0.35,  0.12) 

Gaussian 

(0.37,  0.42) 

(0.31,  0.28) 

(0.3,  0.24) 

Parabolic 

(0.38,  0.22) 

(0.35,  0.15) 

(0.32,  0.18) 

Bessel 

(0.38,  0.2) 

(0.34,  0.19) 

(0.32,  0.18) 
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Inserting  approximation  (22a)  into  emittance  definition  (20a)  produces 


Er  =  a  - ~ - 

fo  3  / 

me  /e 


(l-fcUb  d-fc)Ib 

- o- -  -  «  - 

me  /e  17  kA 


(22b) 


where  a  is  a  dimensionless  coefficient  that  depends  on  the  beam  profile. 
For  example,  for  a  flat-topped  profile, 


-v’’  ■  *  m’  }*■ 


(23a) 


5  (1+g^) 


„  R  (1-f  )IK 

/V  2go  *A 


Jdx 


(23b) 


where  x  =  r  /R  .  Performing  the  integrals  and  substituting  the  results 
into  Eq.  (20a)  yields 

R  (1-fc}Ib  f<1+2gl)  ln<1+2gi>  "  2*l 

efo  =  9l/2  “  3.  U2Z 

2  g  g.  me  /e  L  1 


2g1  -  ln(l+2g1) 


2nl/2 


(23c) 


Using  the  values  for  go  and  g^  from  Table  1,  we  compute  the  emittance 
coefficient  for  a  flat-topped  beam  to  be 


a  -  0.1, 


(23d) 


with  a  weak  dependence  on  b/r^. 

Similarly,  for  the  Gaussian  profile, 


((1-f  )I  )2  22  -1  (k/R) 

«r/V2>  =  (x-e-^*2)  Jdx 

l  goiA  )  K  0  <1+gix) 


(24a) 
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(b/R) 


o  .  R  d-fc)Ib  (  _b2/R2rl  V"' 

<r  At>  ■  k  (1_  c  )  lix 


X  -X 
l+g1X  e 


(24b) 


For  b  >  3R,  Eqs.  (20a),  (24a)  and  (24b)  can  be  combined  into 


R 


'fo  "  Sns 


o6l  me  /e 


1/2 

(gj_  -  d+gj)  y(g2)  -  [y(gx)]2)  ,  (240 


where 


y(g)  =  e 


=  U  £ 


[dt 

td 


can  be  expressed  in  terms  of  an  exponential  integral.  Typically,  the 
emittance  coefficient  for  a  Gaussian  beam  is  given  by 


a  =  0.5, 


(24d) 


five  times  larger  than  the  flat-topped  coefficient  (23d). 
For  the  parabolic  profile, 


«'/V2>  -  !  ^ 

l  Va  J  J  (i+gj_ 


=*m 
... 2  » 


(25a) 


„  2  /jC  v  2R  (1  fc^Ib 
<r  ^  ‘  “ IA~ 


rdx  jso^3> 

J  1+giX 

r\  J- 


(25b) 


Performing  the  integrals  and  substituting  the  results  into  Eq.  (20a)  yields 
1/2  R  (l-fc)Ib  r 

efo  =  I  -372  -37^  <2+3gi>  ln<1+3gi>  -  6gi 
gogl  mc  /e 


2,1/2 


-  I  gl3  (I  gl  +  3gl  ~  (1+3gl}  ln<1+3gi>)  ]•  <25c> 


A  typical  emittance  coefficient  for  a  parabolic  profile  is 
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a  =  0.2,  (25d) 

roughly  half-way  between  the  flat-topped  and  Gaussian  coefficients.  A 
Bessel  profile  should  produce  similar  results. 


C.  Consequences  of  Emittance  Growth 

The  preceding  equations  demonstrate  that  the  emittance  increase  e^o 
from  foil  focusing  scales  linearly  with  beam  current  and  radius,  is 
independent  of  y,  and  depends  strongly  on  the  beam  profile.  High-current 
beams  with  broad  radial  wings,  such  as  Gaussian  or  Bennett  profiles,  suffer 
the  largest  increase.  As  Eqs.  (23d)  and  (24d)  illustrate,  the  presumption 
of  a  flat-topped  profile  can  underestimate  e^o  by  a  factor  of  five  or  more. 

The  significance  of  can  be  determined  by  comparing  it  with  the 
emittance  desired  for  the  beam  outside  the  accelerator.  A  paraxial  beam 
propagates  in  equilibrium  in  a  gas  only  if  its  emittance  equals 

eeq  ■  VR  <Ic££"A>1/2-  <26> 


where  <  1^  is  an  effective  current  that  characterizes  the  pinch  force 

felt  by  the  beam.  In  a  dense  gas,  I  ff  equals  the  sum  of  the  beam  current 

err  i  q 

and  the  inductively  generated  plasma-return  current.  In  a  low-pressure 

gas,  I  ^  measures  the  strength  of  the  plasma  ion  channel.1^ 

The  cumulative  emittance  gained  by  a  beam  after  passing  through  n 

foils,  each  with  an  emittance  increase  of  tc  ,  is  less  than  e  only  if 

fo’  eq  J 

r.  <  - - f  A  ?)  (27] 

l«(l-fc)Ib]2 


assuming  constant  radius  R.  For  f  =  0  and  I  ~  1^,  the  number  of  foils 
is  limited  to  n  <  I^/a^I^,  which  depends  strongly  on  beam  profile  through 
the  parameter  a.  To  increase  the  number  of  allowed  foils,  or  to  transport 
high-quality  beams  with  1^^  «  1^,  plasma  must  be  added  to  produce  f  >  0. 
Adding  plasma  eliminates,  however,  the  need  for  foil  focusing.  Emittance 
growth  from  anharmonic  focusing  thus  limits  the  utility  of  foil  transport, 
especially  for  high-current  beams  with  broad  radial  wings. 
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The  emittance  increase  e^o  could,  if  desired,  be  used  to  alter  the 
emittance  envelope,  e(C),  of  an  initially  cold  beam.  For  example,  a  beam 
of  constant  radius  but  rising  current  would  suffer  the  greatest  emittance 
increase  in  the  beam  body,  regardless  of  variations  in  y*  Alternatively, 
e(C)  could  be  modified  by  varying  the  bean  radius  R  before  the  beam  strikes 
the  foil. 

We  have  thus  far  neglected  changes  in  the  beam  profile  at  the  foil. 

Such  changes  either  increase  or  decrease  the  emittance  through  the 
parameter  e ^  defined  in  Eq.  (20b).  Because  an  anharmonic  lens  itself 
causes  the  beam  profile  to  vary,  is  often  comparable  in  magnitude  to 
E£q.  The  cumulative  effect  of  from  many  lenses  is  usually  small, 

however,  because  of  the  varying  sign  and  lack  of  coherence  of  from  one 
lens  to  another;  see  Sec.  IV-C.  Not  included  in  our  analysis  is  the 
emittance  change  produced  by  the  anharmonic  self-pinch  force,  proportional 
to  (1-f  )I,  ,  between  foils.  In  vacuum,  the  repulsive  self-force, 

CD  2 

proportional  to  I^/y,  alters  the  emittance. 

IV.  MULTI-FOIL  TRANSPORT 


A.  Closely-Spaced  Foils 

The  previous  sections  dealt  with  focusing  from  a  single  foil.  A  series 
of  foils  raises  new  issues,  one  of  which  is  interactions  between  adjacent 
foils.  Expansion  (3)  shows  that  foils  interact  electrically  when  they  are 
separated  by  less  than  the  pipe  diameter,  2b.  To  evaluate  these 
interactions,  consider  a  foil  located  at  z  =  0  and  an  adjacent  foil  located 
at  z  =  d.  The  foil  potential  for  0  <  z  <  d  can  then  be  expressed  as 


$f(r,z)  =  £  Jo(Xnr/b)  [Ar  coshO^z/b)  +  Bn  sinhO^z/b)  J ,  (28) 

n=l 

in  place  of  Eq.  (3).  Here  Ar  is  given  by  Eq.  (4),  and  8n  is  given  by 


2  csch(Xd/b)  ur  n 

3_  =  -  - 0 -  Jdr  r  JQ(Xor/b)  [± 


n 


*n2jl<V  0 


—  r  E 
3r  ro 


(r,d)’ 


-  Ar  coth(Xnd/b). 


(29) 
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Note  that  the  expansion  coefficients  Ar  and  Bn  are  determined  solely  by  the 

ambi'ent  field  E  on  each  of  the  two  foils.  A  similar  but  more  restrictive 

ro  .10 

solution,  involving  a  double  summation,  was  given  by  Humphries. 

If  the  beam  changes  little  between  the  two  foils,  the  foil  cell  behaves 

as  a  single  thin  lens  with  Ero(r,d)  =  ErQ(r,0)  The  momentum  increase 

imparted  within  the  foil  cell  is  then  given  by 


&Pr(r)  = 


-  jdz 


c  3r 


An  J^r/b) 


coshCXjjd/b)-! 

sinhC^d/b)  * 


(30) 


with  Ar  given  by  Eq.  (7).  A  series  of  foil  cells  thus  behaves  as  a  series 
of  thin  lenses,  each  producing  a  deflection  angle  that  takes  the  same 
form  as  that  for  an  isolated  foil,  except  that  the  terms  in  the  sums  over  n 
must  be  multiplied  by  [cosh(XI)d/b)-l]/sinh(Xnd/b) .  The  thin-lens 
(paraxial)  treatment  is  valid  in  this  case  only  if 


ft  »  d,  ■  (31) 

where  f^  is  to  be  computed  from  definition  (12)  and  result  (30).  This 
condition  replaces  the  previous  paraxial  condition,  f^  »  b,  for  widely 
spaced  foils. 


B.  Axisymmetric  Stability 

An  important  issue  in  multi-foil  transport  is  beam  stability.  To  study 
axisymmetric  stability,  we  treat  the  foils  as  ideal  lenses,  3f^/3r  =  0,  and 
assume  that  the  beam  propagates  force-free  in  vacuum  (f£  =  0)  between 
foils.  The  dependence  of  f^  on  beam  radius  R  is,  however,  retained. 

For  thin  lenses,  the  particle  orbits  are  given  by 

x  -  =  x  +  x'  d,  (32a) 

n+1  n  n 


and 


x'  =  x' 


-  *A 


(32b) 


Here  x  is  the  particle  position  at  the  nth  lens,  x'  =  v  /c  is  the  particle 
n  n  x 
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velocity  just  after  the  nth  lens,  and  f  is  the  focal  length  of  the  nth 

len^.  Like  equations  describe  motion  in  the  y-direction. 

For  constant  lens  spacing  d,  we  can  eliminate  x^  to  obtain  a  third- 

order  difference  equation  for  x  : 

n 


x  -  -  2(1  -  d/2f  )x  +  x  -  =  0. 
n+1  •  '  n'  n  n-1 

For  constant  fn,  this  equation  predicts  unbounded  growth  unless 


(33) 


d  <  4fn.  (34) 

Exponential  growth^  occurs  at  d  >  4f  ,  while  linear  growth  occurs  at 
d  -  4fn. 

For  foils,  the  focal  length  f  is  not  constant  but  varies  linearly  with 
the  rms  beam  radius  Rn«  This  dependence  appears  to  be  stabilizing,  because 
condition  (34)  is  ultimately  satisfied  as  Rr  becomes  large.  More  careful 
analysis  reveals,  however,  that  the  dependence  on  Rn  is  destabilizing  until 
Rn  becomes  very  large.  To  show  this,  we  replace  difference  equation  (33) 
with  an  equivalent,  fourth-order  expression  for  R  as  derived  in  the 
appendix: 


<1-d/2£n)R„'2 

+  [l-4(1-d/2£n+l><1-d/2fn>lI<1-d/2f„+1)Rn?i  '  <1-d/2£n>R„2l 

-  (l-d/2fn+1)Rn^1  =  0.  (3b) 


Stability  requirement  (34)  for  constant  f  can  be  confirmed  from  this 

expression  using  solutions  of  the  form  R  =  R  exp(iwn)  and  R  =  nR  . 

no  no 

For  f  «  R  ,  Eq.  (35)  is  nonlinear  and  difficult  to  analyze.  We 

therefore  consider  a  small  perturbation  about  an  equilibrium  solution  of 

constant  radius  R  : 

o 


Rn  =  Rq  [1  +  S  exp(iwn) ] , 


(36) 


with  &  «  1.  Inserting  this  form  into  Eq.  (35)  produces  to  order  8, 
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2(l-k)  exp[ico(n+2)]  +  k  exp(iwn) 

+  ( 2— k) [  1  -  4(l-k)2]  (exp[i«(n+l) ]  -  exp(icon)} 

-  2(l-k)  exp[iw(n-l)]  -  k  exp[iw(n+l) ]  =  0,  (37a) 

where  k  =  d/2f  is  a  geometrical  factor  and  f  =  f  (R  )  is  the  nominal 
'o  o  n  o  i 

focal  length. 

Equation  (37a)  can  be  reduced  to 

cos2(co/2)  =  (1— k)  (1— k/2) ,  (37b) 

which  indicates  that  stable  solutions  (real  w)  are  possible  only  if  k  <  1. 

One  can  similarly  show  that  solutions  of  the  form  R  =  R  [1  +  (— )nn6]  are 

no 

unstable  for  k  =  1.  We  thus  conclude  that  stability  is  possible  in  the 
linearized  regime  only  if 

k  =  d/2f  <  1.  (38) 

o 

This  condition,  which  applies  when  f  <=  R  ,  is  twice  as  restrictive  as 

n  n 

condition  (34)  for  constant  fR. 

If  condition  (38)  is  not  satisfied  at  injection,  Rr  and  f  grow  until 
the  condition  is  satisfied.  This  suggests  that  unstable  solutions  would 
saturate  if  the  beam  were  to  expand  slowly  from  foil  to  foil.  To  test 
whether  condition  (38)  is  actually  necessary,  we  solved  difference  equation 
(35)  numerically,  using  a  10%  variation  in  Rn-  Saturation  occurred,  if  at 
all,  only  after  the  beam  expanded  to  a  radius  much  greater  than  Rq.  We 
therefore  conclude  that  multi-foil  transport  is  in  general  stable  only  if 
condition  (38)  is  satisfied  at  beam  injection. 

There  are  several  caveats  to  the  stability  analysis.  First,  the 
underlying  assumption  of  paraxial  foil  focusing  is  often  poorly  satisfied 
when  stability  condition  (38)  is  met  only  marginally;  preliminary  numerical 
simulations  indicate,  however,  that  this  criterion  remains  roughly  intact 
even  when  the  paraxial  approximation  fails.  Second,  the  analysis  neglects 
emittance  growth  from  anharmonic  foil  focusing  and  scattering.  Emittance 


16 


growth  proportional  to  radius  Rn  ultimately  causes  "stable"  beams  to  expand 
lindarly  with  foil  number  n;  to  minimise  n,  the  foil  spacing  d  should  be 
maximized,  subject  to  stability  condition  (38).  Third,  the  analysis 
assumes  constant  beam  energy  y  (and  foil  spacing  d).  If  r  were  to  double 
over  ten  or  fewer  foils,  beam  stability  and  e  ittance  growth  become  much 
less  troublesome;  practical  limits  on  average  accelerating  gradients  make 
this  requirement  difficult  to  achieve,  however,  for  high-current  beams. 

Pnd  fourth,  the  analysis  ignores  self  forces,  plasma  forces,  and  current 
loss  to  the  pipe.  In  practice,  instability  often  manifests  itself  not  as 
total  disruption  but  as  gradual  loss  of  the  beam  to  the  pipe  walls. 

C.  Matching 

Previous  investigators^’ ^  chose  foil  spacing  based  on  <\  matching 

criterion  in  which  the  beam  propagates  unchanged  from  one  fx.il  cell  to  the 

next.  Here  we  derive  a  matching  criterion  and  determine  if  it  is 

compatible  with  stability  condition  (38). 

For  a  matched  beam,  the  minimum  radius,  R  .  ,  occurs  at  the  mid-point 

mm 

of  each  foil  cell,  while  the  maximum  radius,  R  ,  occurs  at  the  ends.  The 

max 

beam  expansion  rate,  R'  =  dR/dz,  is  thus  anti-symmetric  about  the  cell  mid¬ 
point  and  about  each  foil.  Each  foil  must  therefore  change  R'  to  -R' : 


S6f  “  -  EmaX/£t  -  -  2R'<a>- 


(39) 


where  R’(d)  is  the  expansion  rate  as  the  beam  approaches  a  foil. 

The  expansion  rate  R'  depends  on  the  forces  acting  inside  the  foil 

cell.  Let  us  consider  the  simple  example  of  emittance-dominated,  free 

12 

expansion  in  vacuum  (f£  =  0).  For  constant  y,  the  envelope  equation 
dictates  that  the  beam  expands  past  the  cell  mid-point,  z  =  d/2,  as 


2  2 
R  (z)  =  R  f  + 

mm 


(  t 


n2 


yR  . 
mm; 


(z  -  d/2)  . 


(40) 


Differentiation  yields  the  expansion  rate  just  prior  to  a  foil: 


R'(d)  = 


yR, 


min 


2R 


(41) 


max 
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Combining  this  with  focusing  requirement  (39)  and  Eq.  (AO)  produces 

i 

n2 


£cd  ■ 


^max^min 


2p  4  «2 

Y  Rmin  ,  d 
=  9  + 

A 

Solving  for  the  foil  spacing  d  yields  the  matching  requirement 


(A2) 


d  =  2f . 


/ 

<1  + 

1  - 

r  ?  1 
Y*min 

2  , 

1/2' 

E  f£  1 

\ 

\  / 

/ 

>  2ft. 


(A3) 


Matching  requirement  (A3)  violates  stability  condition  (38),  indicating 
that  stable,  matched  propagation  is  not  possible  with  foils.  That  is,  beam 
stability  demands  closely  spaced  foils  that  overpinch  the  beam  and  cause  it 
to  vary  from  cell  to  cell.  Fluctuations  in  beam  radius  need  not  disrupt 
the  beam,  but  they  do  make  controlling  it  more  difficult  and  may  accelerate 
emittance  growth.  These  fluctuations  also  reduce  contributions  from  the 
emittance  change  by  destroying  beam  coherence  from  cell  to  cell.  Note 
that  stable,  matched  propagation  is  possible  for  thin  lenses  of  constant 
focal  length,  according  to  stability  criterion  (3A). 


V.  OTHER  CONSIDERATIONS 


A.  Foil  Scattering 

Beam  electrons  scatter  as  they  pass  through  a  foil.  Although  energy 

loss  from  inelastic  scattering  is  usually  negligible  for  the  beam,  it  heats 

and  may  melt  the  foil.  Foil  durability  is  determined  by  the  thermal 

properties  and  energy  loss  function  of  the  foil. 

Elastic  scattering  is  often  a  greater  concern  and  generates  an  angular 

spread  in  the  beam  of  characteristic  width  50g  a  Z/t^/y,  where  t r  is  the 

foil  thickness  and  Z  is  the  atomic  number  of  the  foil  material.  This 

2  2 

angular  spread  increases  the  emittance  squared  by  an  amount  e  =  (50  yR)  , 

s  s 

independent  of  beam  current  and  energy.  Unless  the  foil  is  a  highly 
transparent  mesh,*^  eg  typically  exceeds  the  increase  from  anharmonic 
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focusing  at  beam  currents  below  17  kA.  Note  that  the  total  emittance 
increase,  +  eg  ,  is  independent  of  beam  energy  y;  an  energy  ramp  does 

not,  therefore,  directly  alter  the  emittance  envelope,  e(£) . 

B.  Focusing  by  Apertures 

An  apparent  means  of  avoiding  or  reducing  collisional  effects  is  to  use 
conducting  apertures.  Apertures  short  much  like  foils,  but  they  are 
ineffective  at  focusing.  To  show  this,  we  use  Gauss's  law  to  express  the 
radial  impulse  delivered  by  an  isolated  foil  or  aperture  as 

CD 

Spr(r)  =  -  |  Jdz  Erf(r,z)  =  -  Jdr'  r'  of(r'),  (44) 

-to  0 

where  a ^  is  the  surface  charge  density  on  the  foil  or  aperture.  In 

deriving  Eq.  (44)  we  have  used  the  fact  that  the  axial  foil  field  E^  =  0 

at  z  =  ±®.  A  net  impulse  6pr  is  thus  delivered  only  to  beam  electrons  that 

encounter  enclosed  surface  charge.  For  beam  electrons  passing  through  a 

hole,  the  enclosed  surface  charge  is  zero.  Hence,  apertures  focus  only 

18 

that  portion  of  the  beam  not  passing  through  the  hole.  Thin  apertures 

might,  however,  be  useful  in  that  they  can  focus  the  expanding  beam  wings 

while  leaving  the  cold  laminar  core  of  the  beam  unperturbed. 

C.  Foil  Focusing  in  a  Plasma 

In  the  preceding  analysis,  we  allowed  for  a  constant,  self-similar 

charge-neutralization  fraction,  fc  <  1.  This  neutralization  reduces  the 

focusing  action  of  the  foil,  and  additionally  keeps  the  beam  pinched 

between  foils.  An  important  example  is  the  ion-focused  regime  (IFR)  in 

which  the  beam  ionizes  a  low-pressure  gas  and  then  electrostatically  expels 

the  plasma  electrons,  leaving  th'  heavy  ions  to  pinch  and  guide  the  beam.^ 

2 

The  simplicity  and  utility  of  IFk  transport  warrants  further  discussion, 
particularly  at  the  transition  from  vacuum  to  gas  where  foils  are  commonly 
employed . 

Near  the  foil,  the  plasma  electrons  are  not  ejected  radially  because 
the  foil  shorts  the  transverse  electric  fields.  However,  the  axial  foil 
field  Ez£  is  strong  in  this  region  and  removes  the  plasma  electrons 
longitudinally,  either  by  pulling  them  into  the  foil  or  pushing  them  into 
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regions  where  they  are  expelled  radially.  In  either  case,  only  the  ions 
remain,  and  hence  the  assumption  of  uniform  f  over  a  foil  focal  length  is 
justified. 

Foils  are  commonly  used  to  separate  regions  of  different  gas  pressure 
and,  consequently,  different  neutralization  fraction  f£.  Because  a 
perfectly  conducting  foil  disrupts  electromagnetic  communication,  the 
fields  and  impulse  can  le  computed  independently  on  each  side  of  the  foil. 
In  the  thin-lens  limit,  the  total  impulse  and  focusing  effect  are  given  by 
averaging  f  on  the  two  sides.  For  example,  at  a  transition  from  vacuum  to 
gas,  the  foil  focal  length  f^  and  emittance  increase  can  be  computed  as 
before,  using  f  set  to  half  its  value  in  the  gas  region.  The  analysis  of 
Sec.  IV-A  for  closely  spaced  foils  should,  however,  be  applied  without 
modification.  Conservation  of  beam  and  field  energy  dictates  that  a  change 
in  f  across  a  foil  alters  the  axial  beam  momentum  and  energy  by  an  amount 
proportional  to  the  change. 


D.  Non-Axisymmetric.  Effects 

In  the  previous  analysis,  the  currents  and  charges  induced  in  the  pipe 
walls  were  axisymmetric  and  had  no  effect  on  the  beam.  Off-axis  beams, 


however,  induce  non-axisymmetric  wall  currents  and  charges  that  deflect  the 

beam.  The  wall  return  currents  magnetically  repel  and  center  the  beam, 

while  the  wall  charges  electrostatically  attract  the  beam.  Because  foils 

neutralize  the  electrostatic  fields,  they  provide  a  centering  impulse,  as 
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well  as  a  focusing  impulse,  to  the  beam.  ’ 

For  off-axis  beams,  the  foil  potential  can  be  expanded  in  terms  of 
trignometric,  exponential,  and  Bessel  functions.  The  centering  impulse  can 
then  be  calculated  using  the  technique  given  earlier  for  the  focusing 


impulse.  For  small  beam  offsets,  y^  «  b,  a  simpler  estimate  is  possible. 
In  this  case,  an  average  foil  centering  force  ~  eCl-f^I^y^/b^c  persists 
over  a  distance  -  2b,  producing"*  a  centering  impulse  Sp^  =  2e(l-fc)I^y|3/b. 
Using  analysis  analogous  to  Eqs.  (32)-(34),  one  can  show  that  the  impulse 


is  stabilizing  only  for  foil  spacings 


d  <  4ybpz/sPy  «  2bIA/(1-fc)Ib* 


(45) 


This  criterion  for  non-axisymmetric  stability  is  generally  well  satisfied 


when  the  axisymmetric  criterion  (38)  is  satisfied  and  b  >  2R. 

5  6 

The  foil  centering  impulse  has  been  claimed  ’  to  help  stabilize  the 

beam.  This  prediction,  like  condition  (45),  is  based  on  a  rigid-rod 

treatment  of  the  beam.  However,  the  finite  transverse  electromagnetic 

propagation  time  (~  2b/c)  of  the  wall  forces  can  cause  the  foil  centering 

impulse  to  become  out  of  phase  with  high-frequency  oscillations  of  the  beam 

centroid.  The  impulse  then  causes  the  oscillations  to  grow  from  one  foil 

to  the  next,  or  even  with  distance  C  into  the  beam  pulse-  Hence,  although 

closely  spaced  foils  damp  low-frequency  transverse  oscillations,  they  may 

produce  instability  at  high  frequencies. 

A  related  concern  is  the  beam  breakup  instability  (BB>')  which  develops 

as  transverse  beam  motion  induces  and  couples  to  electromagnetic  waves 
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within  the  foil  cells.  ’  *  Adding  foils  increases  the  number  of  foil 
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cavities  and  likely  enhances  BBU  growth.  Colombant  and  Lau  have  recently 

concluded  that  BBU  could  explain  the  poor  transport  observed  at  high 

q 

current  in  the  multi-foil  experiments  of  Meger.  Poor  transport  is  also 
predicted,  however,  by  the  axisymmetric  analysis. 

Higher-order  azimuthal  magnetic  instabilities  may  develop  as  the  foils 
partially  neutralize  the  electrostatic  fields.  Humphries  and  Ekdahl^ 
concluded  from  foil  experiments  in  vacuum  that  localized  magnetic  pinching 
causes  beam  filamentation.  This  azimuthal  bunching  saturates  by  raising 
the  transverse  beam  pressure  and  emittance.  Emittance  growth  was 
particularly  pronounced  for  initially  hollow  beams. 

VI.  CONCLUSION 

5  10 

The  theoretical  analysis  presented  here  supports  earlier  work  ’ 
indicating  the  utility  of  foil  focusing  of  electron  beams,  but  it  also 
indicates  potential  pitfalls  and  limitations.  Problems  arise  because  foil 
focusing  depends  on  the  beam  parameters  and  because  it  is  anharmonic.  For 
example,  the  dependence  of  the  foil  focal  length  on  the  beam  current, 
radius,  and  energy  not  only  makes  foil  focusing  difficult  to  control,  as 
different  beam  segments  respond  differently,  but  it  also  leads  to  the 
surprising  conclusion  that  matched,  stable  propagation  is  not  possible. 
Instead,  the  beam  radius  oscillates  in  an  unpredictable  although  mild 
manner  as  the  beam  propagates  from  foil  to  foil.  These  oscillations 
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aggravate  emittance  degradation  from  anharmonic  foil  focusing  and 
scattering.  For  beam  currents  comparable  to  the  Alfven  current,  emittance 
degradation  from  anharmonic  focusing  is  likely  to  become  excessive  after 
ten  or  fewer  foils.  In  general,  we  find  that  multi-foil  transport,  like 
other  techniques,  works  best  for  beam  currents  that  are  small  compared  with 
the  Alfven  current. 
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APPENDIX 


The  transport  properties  of  a  series  of  thin  lenses  at  fixed  spacing  d 
are  described  by  Eq.  (33).  Here  we  derive  an  equivalent  expression  for 
ideal  lenses  in  terms  of  the  rms  beam  radius,  Rn>  The  latter  expression  is 
useful  when  the  focal  length  f  varies  with  R  . 

Because  Rn  involves  an  average  of  x^  plus  y ^ ,  where  (x  ,y  )  is  a 
particle  position,  we  begin  by  squaring  Eq.  (32a)  for  x  ,  to  obtain 


2x  x'd 
n  n 


+  (x^d)-. 


(Al) 


We  next  multiply  Eq.  (32a)  for  xr+1  times  Eq.  (32b)  for  x^+1: 

9  x  ,  (x„  +  x'd) 

Vrf.l  ■  Vn  *  <■*'«>  *  -  -B±SS - “ 


VA  +  <*A)2d 


=  [l  -  - 1  X  X'  - 

*■  w  "  ” 


n+1 

(V  XAd)xn  *  W*A.l*  VA.I  >d 

£„.l 


fn.l  n+1  n+1  ‘  £n+l 


n+1 


l^n+1 


+  (x£)M. 


(A2) 


Here  we  have  used  Eqs.  (32a)  and  (32b)  to  eliminate,  in  selected  places, 
xn+l  and  xn*  Eliminating  x^2  from  Eqs.  (Al)  and  (A2)  produces 


(x 


n+1  n+1 


x  x/ )dl 
n  n 


(A3) 


We  now  square  Eq.  (33),  written  in  the  form  xfl  2  =  2(l-d/2fn)xn 


x 


n+15 


(A4) 


where  we  have  again  used  Eq.  (32a)  to  eliminate  xr+1  from  the  last  term. 
After  rearrangement,  this  equation  reduces  to 
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(A5) 


x  x' 
n  n 


xn+l  ~  xn-l 
4(1  -  d/2fn)d 


Result  (A5)  can  be  used  to  eliminate  xn+ix„+i  an<^  xnxn  from  Eq.  (A3), 
producing 


*  [l-4(l-d/2£n+1)(l-'i/2fn)]I(l-d/2£ntl)xn21  -  (l-d/2£n)xn2] 

-  (l-d/2fn+1)xn^1  =  0.  (A6) 

A  similar  expression  applies  to  y  .  If  we  sum  and  average  these 
expressions,  and  invoke  the  perfect-lens  assumption  that  f  is  independent 
of  the  particle  positions  (xn>yn) »  we  obtain  Eq.  (35): 

<1-d/2£„>En!2 

+  [l-4(l-d/2fn+1)(l-d/2fn)][(l-d/2fn+1)RnJ1  -  <l-d/2fn)Rn2] 

-  (l-d/2fn+1)Rn21  =  0,  (A7) 

2  2  2 
where  R  is  the  average  of  x  +  y  . 
n  °  n  vn 
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Fig.  1  The  foil  focal  length  f^  expressed  in  units  of  Rl^^l-f^JI^.  No 
beam  electrons  exist  outside  r/R  =  /2  for  the  flat-topped  profile  (la), 
outside  r/R  =  /3  for  the  parabolic  profile  (lc),  and  outside  r/R  =  1.8  for 
the  Bessel  profile  (Id).  The  Gaussian  profile  (lb)  extends  to  the  wall. 
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